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Abstract 

We study the relationship of the spectral form factor with quantum as well as classical probabili- 
ties to return. Defining a quantum return probability in phase space as a trace over the propagator 
of the Wigner function allows us to identify and resolve manifolds in phase space that contribute 
to the form factor. They can be associated to classical invariant manifolds such as periodic orbits, 
but also to non-classical structures like sets of midpoints between periodic points. By contrast to 
scars in wave functions, these features are not subject to the uncertainty relation and therefore 
need not show any smearing. They constitute important exceptions from a continuous convergence 
in the classical limit of the Wigner towards the Liouville propagator. We support our theory with 
numerical results for the quantum cat map and the harmonically driven quartic oscillator. 
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Introduction Evidence abounds that the spectrum of quantum systems bears information 
on the corresponding classical dynamics, in particular on manifolds invariant under time evo- 
lution. The Gutzwiller trace formula l|] and its numerous ramifications feature specifically 
the set of isolated unstable periodic orbits of classically chaotic systems. The discovery that 
energy eigenfunctions are typically "scarred" along such orbits ^] required to modify the 
picture of ergodic eigenstates and allowed for the first time to directly visualize the impact of 
classical invariant manifolds on quantum mechanical distributions defined on configuration 
or phase space The influence of classical invariant manifolds on time-domain features 
has mainly been studied in the spectral form factor. It inherits its relation to periodic orbits 
from the underlying spectral density via the Gutzwiller trace formula. Being bilinear in 
the spectral density, it involves pairs of orbits and their interfering contributions. A host 
of research work has been dedicated to evaluating the double sum over periodic orbits that 



ensues 



Q. 



Only recently, the full sum could be tamed, thus providing an exact semiclassical 
account of the form factor j^. 

A step towards more global and immediate relationships to the classical dynamics has 
3een made in the context of the spectral analysis of systems with dynamical localization 
6, 7], in the form of a direct relation of the spectral form factor K{t) jy, Q] with the classical 
probability to return P^^^(t). For chaotic systems it reads 

K{t) ^ {2/f3)rP:^,{tnr), (1) 

where /3 = 1 for systems invariant under time reversal and 2 otherwise. Being based on 
the diagonal approximation, the expression is valid for times short compared to the Heisen- 
berg time tu. A similar relation holds for integrable systems, but without the prefactor r. 
Equation ([T]) calls for a deeper understanding and analysis beyond its original application 
and derivation from the Gutzwiller trace formula, to explore its potential as an alternative 
semiclassical route to spectral analysis. 

In this Letter, we study the relation of quantum and classical return probabilities in 
phase space with the spectral form factor in the light of recent progress in semiclassical 
approximations to the Wigner propagator j^, 9]. This approach has the special merit that the 
interference of orbit pairs is already implicit in quantum return probabilities. They can be 
expressed, like their classical analogues, as traces (not traces squared!) over a corresponding 
propagator, resulting in very direct quantum-classical relations on the same footing. 
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Before tracing, the diagonal propagator of the Wigner function, through its exphcit de- 
pendence on phase-space coordinates, allows to resolve the manifolds in phase space behind 
the contributions to the form factor. Expressing it semiclassically in terms of orbit pairs, 
it turns out that besides the classical invariant manifolds also sets of midpoints between 
them contribute. Hence classical and quantum return probabilities generally cannot coin- 
cide. This implies severe restrictions to the convergence of the Wigner propagator towards 
the classical (Liouville) propagator, at least for the diagonal propagator near such midpoint 
manifolds. That these dominant features of the diagonal Wigner propagator, classical as well 
as non-classical ones, occur in a time-dependent distribution function suggests calling them 
"time-domain scars". By contrast to scars in eigenfunctions, they are not affected by the 
uncertainty relation and therefore allow for an unlimited resolution of classical structures. 

Classical and quantum return probabilities In quantum mechanics, a probability to return 
is generally defined like an autocorrelation function: Introduce a return amplitude a^etit) — 
J qo{q{t)\qo) with |q(t)) = C/(i)|qo), U{t) the time-evolution operator, and square. 



By contrast, a classical return probabihty in phase space is constructed as follows: Prepare a 
localized initial distribution Pro{i^: 0) = (5A(r — Tq), (^aIi") a strongly peaked function of width 
A and r = (p, q) a vector in 2/-dimensional phase space. Propagate it over a time t and over- 
lap it with the initial distribution. The resulting pf^^{ro,t) — J d^-^'r pro(r, i)pro(r, 0) can be 
interpreted as a probability density to return. Here, the time-evolved distribution is obtained 
from the Liouville propagator G"^^(r", t; r', 0) as pr(,(r",^) = / d^-^r' G"^^(r", t; r', O)pro(r', 0). 
Tracing over phase space yields the return probability -Pret(^) — I ^^^''^oPrlti^o,t)- Replacing 
the initial distribution hy S{r — Tq), we have 



To avoid divergences in particular at t = 0, the phase-space integration has to be restricted 
to a finite range AE in energy, if it is conserved, by introducing some normalized energy 
distribution p{E). 

In quantum mechanics, the Wigner function allows for a similar construction. Be- 
ing related to the density operator p{t) by an invertible transformation, W{r,t) — 
J q' e"^^'^' {q + q!/2\ p{t) |q — qV2), its propagator is defined as the kernel that evolves 



a,et{t)f ^\trU{t) 



(2) 




(3) 
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it over finite time, W{r", t") = J d^-^r' Gw(r", t"; r', t')W{r', t'). By analogy, we thus arrive at 
a quantum-mechanical quasi-probability density to return in phase space Pr^{ro,t) = 
Gw(ro, t] To, 0), and a return probability 

PrTit)= /d2^roGw(ro,t;ro,0). (4) 

The integration across the energy shell produces a factor D-h = AE/{d), the effective di- 
mension of the Hilbert space Ti, (d) denoting the mean spectral density. 

Equations (jl]) and ([2]) are equivalent, as becomes clear if we express the propagator of 
the Wigner function in terms of the Weyl propagator, U{v,t) = J d-^g' e~^P'^'/^ (q + q'/2| 
Uit) |q-q72), 

Gw(r",t;r',0) = J d^^r e^(^"-'-')^^f/*(r_, t)t/(r+, t), (5) 

with r± = (r' + r" ± r)/2. Substituting in Eq. (jlj) and transforming to r'^ = tq ± r/2, the 
two integrals factorize, PX(i) = J d'^^r'_U*{r'_,t) J d'^h'_^U{r'_^,t) = \tTU{t)\\ 

Form factor and diagonal propagator Also the form factor is related to the trace-squared 
of the time-evolution operator, K{t/t]i) = |trt/(t)p for t ^ tu/D-^, where tu = h{d). 
The factor normalizes limT-^oo -^(t) = 1- By comparison with Eqs. ([2]) and (HI), 

P,rW = / d2^rGw(r,t;r,0) = DnK{t/tn). (6) 

This remarkable relation expresses the form factor as the trace over a quantity with a close 
classical analogue, not as a squared trace. It is an exact identity and does not involve any 
semiclassical approximation. 

Contrast Eq. ([6]) with ([T]). Both relate K{t) with a return probability, but there is a clear 
discrepancy, manifest in the factor r that appears only in ([T]). This may not be surprising 
given that the two relations refer to return probabilities on the quantum and the classical 
level, respectively. However, if we take into account also Eqs. (131) and (HI), we face a dilemma: 

n fi n ^ 

There is ample evidence [8l, l9j, lll[ that the Wigner propagator generally converges in the 
classical limit to the Liouville propagator, 

lim Gw{r", t; r', 0) = G'\r" , t; r', 0). (7) 

h—*o 

For up to quadratic Hamiltonians, is even identical to it. Were Eq. (JTj) correct also for 
r' = r" — and on the diagonal the Wigner propagator should behave more classically than 
elsewhere — then lim/^^o -Pre™(^) = -^reU^) should hold as well! 



The derivation of Eq. ([T]) suggests that the factor r arises as a degeneracy factor due 
to the coherent superposition of contributions from different points along a given periodic 
orbit, each of which can be interpreted as a periodic point of its own, r measuring the 
magnitude of this set in phase space. We therefore suspect that Eq. ([7j) might fail in the 
presence of constructive quantum interference. This can be substantiated taking into account 
semiclassical approximations for G'w(r", t; r', 0) based on pairs of classical trajectories j^, [o] 
r'^(t), r'j'(t), chosen such that for their respective initial points r^, r' = (r'_ + r'_^)/2, and 
likewise for r'^. Specifically for the diagonal propagator, this requires that both r^}(t) and 
r'^(t) be periodic orbits. The set of midpoints f{t) = {r'^}{t) + r'^(t))/2 then forms a closed 
curve in phase space as well and contributes to the diagonal propagator hence the form 
factor, but need not consist of periodic points proper. 

It is tempting to interpret also the prefactor 2//3 in Eq. ([1]) as a degeneracy factor and to 
look for phase-space manifolds that in time-reversal invariant systems contribute the extra 
weight to P^^{t): They can be found in sets of midpoints between symmetry-related pairs of 
periodic orbits, located in the symmetry (hyper)plane p = 0. Similarly, other non-diagonal 
contributions to the form factor sl can be associated to non-classical enhancements of 
the diagonal Wigner propagator. 

Examples In order to render our argument more quantitative, we first discuss the case 
of discrete time: Consider a set of periodic points rj{n + Nj) = rj{n), n = 0, . . . , Nj — 1, 
of a symplectic map Ai. In their vicinity, the semiclassical Wigner propagator is given by 
Gwi(r", Nj] r', 0) = 6{r" - M^r'), Mj denoting M^^ linearized near r', r". Define midpoints 
fj{m,n) = (rj(m) + Yj{n))/2 (cf. Fig. [T]). By construction, fj(m + Nj,n) = fj{m,n), but 
generally A4^^fj{m,n) ^ fj{m,n). For r' ^ r" ^ fj{m,n), the Wigner propagator carries 
an additional oscillatory factor, 

Gwj(r", Nf, r', 0) = 26{r" - M^r') x 

cos {{jjin) - r,(m)) A (r" - Y')/h). (8) 

From here, tracing reduces to equating r' with r" and summing points. There are Nj periodic 
points on the orbit and Nj{Nj — 1) midpoints (fj(m, n) and fj{n,m) count separately), 
resulting in a total return probability 

PX,(iV,) = N]/\det (M, - /)| = N,P:'{N,). (9) 
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^f-^- r(2,3) 

FIG. 1: Schematic drawing of a set of periodic points with period 5 of a symplectic map with their 
midpoints. 
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FIG. 2: Diagonal Wigner propagator G^f/ {r , n; r , 0) for the quantized Arnol'd cat map at n = 1 
(a) and n = 3 (b). Symbols x, + mark periodic points of the corresponding classical map and 
their midpoints, respectively (for better visibility of the data, symbols have been suppressed in the 
upper half of panel (b)). The Hilbert-space dimension is D-^ = 60. Color code ranges from red 
(negative) to blue (positive). 

The midpoints' contribution thus is responsible for the extra factor r, i.e. here, Nj and 
explains the discrepancy between classical and quantum return probabilities. 

As an example, consider the Arnol'd cat map. It is defined on a torus, r" = Tr'(mod 1), 
r G [0, 1)^, T a 2 X 2 matrix with integer coefficients. We choose the simplest combination 
that allows for quantization 12|, T = (2, 1;3,2). The topology of the underlying classical 
space implies that both position and momentum be quantized, leading to a finite Hilbert- 

unction can be adapted to this discrete 



space dimension D-}{. The definition of the Wigner 

periodic Hilbert space to avoid redundancies Q, 1^. In Fig. [21 we show the diagonal 
Wigner propagator after 1 and 3 iterations of the quantum map. The peaks of the diagonal 
propagator coincide perfectly with the periodic points of the classical map. Moreover, they 
appear with almost single-pixel precision. While the uncertainty relation requires a minimum 
area of Dy^ pixels, this is perfectly admissible for the propagator. To check Eq. we 
compared the trace of the diagonal propagator to analytical results for J2j A^j^/|det {Mj — I) \ 
(2.0 and 50.0, resp.), and found coincidence up to 6 digits. 



FIG. 3: Surface formed by midpoints of a fictitious periodic orbit that is not circularly symmetric 
nor confined to a plane in phase space. It exhibits self-intersections but retains the topology of a 
closed two-dimensional ribbon, see text. 



Going to systems in continuous time, a periodic orbit rj{s) = rj{s + Tj) gives rise to 
midpoints rj(s', s") = (rj(s') + rj{s"))/2. This replaces Eq. (IHl) with 

Gwj(r", t; r', 0) = 26{r" - Mjv') x 

cos ((r,(s") - r,(s')) A (r" - r')IK)5{t - Tj). (10) 

The midpoints now merge into a continuous two-dimensional surface Sj parameterized by 
{s',s"), < s',s" < Tj', the length of the orbit. Topologically it forms a closed ribbon. 
As a consequence, the diagonal propagator consists of a 5-function only in the subspace 
orthogonal to Sj, G^j^r, t; r, 0) = 6{rj_)S(t — Tj)/|det {Mj± — /)|, where Mjj_ is the stability 
matrix restricted to the (2/ — 2)-dimensional subspace r^. Upon tracing, the integration 
over Sj yields a factor Tj^, its effective area, 

PrTjit) = AETf 5{t - T,)/27r;i|det (M,^ - /)|. (11) 

In Cartesian phase-space coordinates r, Sj may have a nontrivial geometry. In general, it 



will exhibit a Wigner caustic 15|, an overlap of three leaves near the center of the orbit, 
owing to the fact that a given point in this region may be the midpoint of more than one 
pair of periodic points on the orbit. The phenomenon can well be observed in Fig. HI If the 
periodic orbit is not confined to a plane, this geometric degeneracy will be lifted, resulting 
in folds and self-intersections, illustrated in Fig. [3] for a fictitious periodic orbit. 

A pertinent example is the harmonically driven quartic oscillator H{p,q,t) = p'^/2m — 
moolq^ /A + m?ujQq^ /QAE\, + Sqcos{ujt + 0) [l^, with generally mixed phase space. In the 
diagonal propagator at t = T = 27t/uj (Fig. H]) we identify a number of isolated peaks at 
periodic points of the classical dynamics, elliptic as well as hyperbolic, and their midpoints, 
and an enhancement over a well-defined region, to be interpreted as the Wigner caustic of a 
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FIG. 4: Diagonal Wigner (a) and Liouville (b) propagators G(r, t\ r, 0) for the harmonically driven 
quartic oscillator at t = T = 2tt/lo, with u)q = 1.0, oj = 0.95, (j) = vr/3, S = 0.07, and = 192.0 
(color code as in Fig. [2]). For better orientation, we superimpose a stroboscopic surface of section 
of the same system (panel (b), black). The figure-cxo structure is the Wigner caustic of a period-T 
torus outside the frame shown (grey). Symbols 0, x mark elliptic and hyperbolic periodic points 
of the classical system, resp., and + their midpoints. 

period-T torus outside the frame shown, as confirms the coincidence with the corresponding 
classical feature in Fig. Hb. 

Refinements and perspectives An alternative access to the Wigner propagator near peri- 
odic orbits is Berry's scar function, a semiclassical approximation to the Weyl propagator 



in the energy domain 17|]. It responds to the special situation close to a periodic orbit j by 
using local curvilinear coordinates: energy, time, and remaining phase-space directions r^x 
perpendicular to the orbit. Transformed to the time domain and substituted for the Weyl 
propagator in Eq. ([5]), it leads to a semiclassical approximation for the diagonal Wigner 
propagator. 



Tj/2Tih 



|det (Mjx -/)| 
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5{v^^)5{t-Tj 



(12) 



The primitive period Tj and the determinantal prefactor measure the length and the effective 
cross section, resp., of the "phase-space tube" around the orbit that contributes to the 
diagonal propagator. By contrast to Eq. (ITUl) . the degeneracy factor Tj appears here already 
before tracing: The use of local coordinates condenses the contributions of periodic points 
as well as midpoints onto the orbit. Equation (1121) does not apply outside the orbit j and 
therefore does not allow for indiscriminate tracing over all of phase space. 

The midpoint contribution to Gw{Y,t]r,0) giving rise to marked non-classical features 
is a manifestation of quantum coherence. It measures the quantum return probability for 
Schrodinger-cat states distributed over different points of the same periodic orbit. In the 
presence of incoherent processes, it decays on the dephasing timescale. The Wigner prop- 



agator, operating on the projective Hilbert space, readily permits including this effect [18 1 
and thus to identify exclusively the classical invariant manifolds, unaffected by the uncer- 
tainty relation, as peaks of a purely quantum-mechanical distribution. Phase-space features 
associated to non-diagonal contributions to the form factor will be even more elusive and 
geometrically more involved, but are in principle accessible to numerical study. 

We have provided analytical and numerical evidence that Eq. ([T]) can be interpreted as 
a global relation between quantum and classical return probabilities which can be broken 
down into contributions of invariant phase-space manifolds. They enter with weight factors 
that measure the size of the set contributing coherently, and lead to important exceptions 
;o Eq. (^^|. Analytical evidence based on presently available semiclassical approximations 
91] indicates they are restricted to the diagonal r' = r" (where they are least expected) and 
hence of measure zero. They are qualitatively different for integrable systems: In action- 
angle variables, the size of the degenerate sets is independent of time and therefore does 
not contribute an extra factor t. This in turn reflects the different dimensions and topologies 
of periodic tori vs. isolated unstable periodic orbits, indicating how to generalize this to more 
involved cases like systems with mixed phase space. Merging the different contributions on 
the classical side into more global quantities like the Frobenius-Perron modes 191] remains 
as a challenge for future research. 
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